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We show that in the regime when strong disorder is more relevant than field quantization the 
superfluid-to-Bose-glass criticality of one-dimensional bosons is preceded by the prolonged loga- 
rithmically slow classical- field renormalization flow of the superfluid stiffness at mesoscopic scales. 
With the system compressibility remaining constant, the quantum nature of the system manifests 
itself only in the renormalization of dilute weak links. On the insulating side, the flow ultimately 
reaches a value of the Luttinger parameter at which the instanton-anti-instanton pairs start to 
proliferate, in accordance with the universal quantum scenario. This happens first at astronomical 
system sizes because of the suppressed instanton fugacity. We illustrate our result by first-principles 
simulations. 
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I. INTRODUCTION 

A quarter-century ago, a quantum phase transition 
from a superfluid (SF) to Bose glass (BG), a many-body 
bosonic counterpart of Anderson localization, was pre- 
dicted theoretically, first in one-dimensional (ID) sys- 
tems by Giamarchi and Schulz (GS) 1 , and later in any 
dimension by Fisher et al. 2 . Over the years, substantial 
progress has been achieved in understanding this fasci- 
nating physics such as describing the interplay between 
disorder and commensurability 3 , including the transition 
to the so-called Mott glass (MG) phase (rather than BG) 
in disordered commensurate systems featuring particle- 
hole symmetry 4 , and the theorem of inclusions, which 
rules out a direct superfluid-to-Mott-insulator transition 
and establishes the Griffiths universality class for tran- 
sitions out of a fully gapped state^. Accurate ground- 
state phase diagrams of the disordered bosonic Hubbard 
model in dimensions d = 2, 3 and partly d = 1 have been 
produced by first-principles Monte Carlo simulations^!, 
which were also used for large-scale simulations of the 
universal critical behavioi^l Nowadays, rapid exper- 
imental developments in the field of ultra-cold atoms 
spark renewed excitementP^HEl 

Although many theoretical challenges of disordered 
bosons have been completed successfully, there remains 
one compelling problem in ID, raised already in Ref. [T] 
Does a strong-disorder critical point exist which is quali- 
tatively different from the one found in Ref. 1 ? That one 
was established by means of a perturbative renormaliza- 
tion group (RG) analysis. It was later proved to be a 
generic scenario for finite disorder strength in Ref. [15] 
while a recent two-loop RG study 16 is consistent with 
the latter conclusion. A strong-disorder critical point 
would also imply the existence of more than one Bose 
glass phase. 



What renders ID systems special, is the role played 
by exponentially rare exponentially weak links (RWLs). 
Even for a classical-field counterpart of the bosonic sys- 
tem, weak links can lead to a critical point at which the 
groundstate loses its macroscopically uniform response 
to a phase-twist applied at the boundaries^. On the in- 
sulating side of this classical transition, a macroscopic 
system becomes equivalent to a single Josephson junc- 
tion with a macroscopically small Josephson coupling. 
Inspired by this observation, Altman et al. proposed 
the strong-disorder scenario for quantum systems based 
on a real-space RG treatment 18 . As argued in Ref. [I9j 
the gradual vanishing of the superfluid stiffness down to 
zero (corresponding to classical criticality) would be in- 
consistent for a quantum system since the proliferation 
of the instanton — anti-instanton pairs is guaranteed in 
the thermodynamic limit when the Luttinger parameter 
reaches the critical value K c (this phenomenon defines 
the generic Giamarchi- Schultz critical pointP^). For the 
SF-BG transition, K c = 3/2, while for the SF-MG tran- 
sition, K c = 2, which is the same as for the SF-Mott 
insulator transition. Thus, the critical value of the Lut- 
tinger parameter in the putative strong-disorder scenario 
has to exceed if c ? which forces the strong-disorder crit- 
icality of Kosterlitz-Thouless type to develop prior to the 
generic mechanism of Refs. H1andfT5l 

In this paper, we show that the RG flow of K{L) can 
retain its classical-field behavior down to the K(po) = 
K c value, where the universal phenomenon of instanton 
proliferation inevitably takes place on sufficiently large 
scales. This scenario occurs when if (£o) = ifo ^> 1 at 
the microscopic cutoff scale ^.e., in the most natu- 
ral regime for the putative strong-disorder scenario. The 
quantum renormalization of the RWLs is irrelevant for 
the flow behavior. In this limit, not only the putative 
strong-disorder critical point is ruled out, but also the 
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FIG. 1. Renormalization group flow of the Luttinger pa- 
rameter which starts as classical flow and crosses over to the 
standard Kosterlitz-Thouless flow at much larger scales if the 
former saturates to a value slightly above or below K c = 3/2. 
These large scales may not be available for numerical analysis. 



FIG. 2. Single-stage renormalization group flow in the hy- 
pothetical strong-coupling scenario. Having passed the corre- 
lation length ro corresponding to the condition K(L = ro) = 
K c > 3/2, the flow rapidly brings the system to localization. 



generic instanton-pair induced flow is practically inoper- 
ational because of the extremely low fugacity of pairs of 
mesoscopic size. As a result, a typical system of linear 
size L that should formally be called an insulator, be- 
cause K(L) < K C1 retains a superfluid response. This 
behavior (illustrated in Fig. [T]) — as we will later see, this 
is the hallmark of our scenario, distinguishing it from 
both the Giamarchi-Schulz case and the quantum strong- 
disorder regime advocated in Ref. [18] — persists over a 
broad range of scales till K(L) « 1 (see further on in 
the text). For the same reason of extremely low instan- 
ton fugacity the critical RG flow features a pronounced 
two-stage character (see Fig. [I] and note that the acces- 
sible sizes may be too small to observe both stages of the 
flow), as opposed to the single-stage flow of the hypo- 
thetical quantum strong-disorder regime (see Fig.|2|. 

We also discuss that a strong-disorder scenario with 
the critical value of K c significantly above the GS value 
may not look quantitatively similar to Fig. [T] because 
if the classical flow levels off at K significantly above 
the 3/2 value, the proliferation of instantons with small 
fugacity (or backscattering processes) becomes impossi- 
ble. This implies that any strong-coupling scenario with 
K > 3/2 should have a single-stage RG flow illustrated 
m Fig. [2] We also prove that any quantum critical flow 
possesses self- averaging properties and thus is subject to 
asymptotic hydrodynamic description. As far as we see, 
this cannot be reconciled with any strong-coupling sce- 
nario. 

While the condition Kq ^> 1 is necessary to guaran- 
tee the applicability of our theory, in practice, the meso- 
scopic classical-field strong-disorder flow can take place 
already at Kq ~ 3, which we clearly see in first-principles 
simulations of a J-current type model up to system sizes 
L ~ 1000. 



The rest of the paper is organized as follows. In Sec. Ill] 



we develop our theory. In Sec. Ill we address the sub- 
tle issues of self- averaging under the conditions of strong 
disorder, arguing that the observables of choice are the 
median values and appropriately defined widths rather 
than the standard cumulants — the expectation value and 
variance, or the dispersion which is the square root of the 
variance. In Sec.[lV|we present the model and the results 
of an illustrative simulation. Our conclusions are formu- 
lated in Sec. |V| 



II. THEORETICAL ANALYSIS 

The main assumption — to be justified later — of the in- 
stanton approach 15 is that up to the critical point the 
system remains a well-defined Luttinger liquid charac- 
terized by a finite parameter K. This allows one to 
describe the coarse-grained behavior of the system with 
Popov's hydrodynamic action where = &(x,y) 

is the (1 + l)-dimensional field of superfluid phase (x 
is the spatial coordinate and y = cr is the imaginary 
time re-scaled with the sound velocity c). The instan- 
tons are the virtual quantum events which result in the 
phase slips of <£>. In (1 + 1) dimensions they are vortices 
in the field $>(x,y). An important distinction between 
the instanton and its 2D vortex counterpart is that the 
contribution of an instanton to the partition-function in- 
tegrand comes with the x-dependent factor e ±ll ^ x \ in 
which the sign is defined by the sense of the vortex ro- 
tation and j(x) — 2tt J n(x') dx' , where n(x) is the 
disorder-dependent microscopic expectation value of the 
number density at point x. In view of the random phase 
factor j(x), the net contribution of instantons to the par- 
tition function survives after integration over x only if in- 
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stantons come in instanton-anti-instanton pairs with two 
x-coordinates being microscopically close to each other 
to ensure that their phase factors compensate each other. 
[This is not the case for the MG where exact particle-hole 
symmetry and commensurability render n{x) = integer 
and the phase factor irrelevant. It is precisely this aspect 
that leads to the difference between the K c values for BG 
and MG transitions.] For our purposes, we do not need to 
consider the general case of RG flow of K(L) due to the 
instanton-anti-instanton (or vortex-antivortex, for brief- 
ness) pairs. When the gas of pairs is extremely dilute, 
the criterion for their proliferation is derived based on the 
statistics of a single pair with a hierarchy of length scales 
starting from the microscopic cutoff £o to the system size 
L being absorbed into the scale dependent value of K(r). 
The proliferation takes place when the "effective energy" 
of the pair (in effective temperature units, T = 1) 

f L Kir) 

£pair = 2 / — ^ dr + E c , (1) 
r 

with E c the (unknown) vortex pair core energy, becomes 
lower than the entropy given by its three degrees of free- 
dom: two in the ^/-direction and one in the x direction 
(there are 4 degrees of freedom for the MG case) 

S B G=31nL, Smg — 4 In L, (2) 

so that the effective free energy turns negative. In the 
genuine thermodynamic limit of L — oo, this happens 
at the universal value K(oo) = K c . For i^(£o) ^ 1 
and given the slow decrease of K with distance, an enor- 
mously large system size is required for the entropy term 
to compete with the contribution in i^ pa ir coming from 
distances where K(r) ^> 1. This explains the extremely 
low fugacity of the vortex-antivortex pairs and accounts 
for the fact that E pa i r (L) can remain larger than S(L) 
even when K(L) is already substantially lower than K c . 

To force a system with large Kq ending up in the in- 
sulating state one needs strong disorder which produces 
enough weak links to renormalize K down to K c . For 
the flow to have classical- field character, the microscopic 
physics of RWL does not necessarily have to be classi- 
cal. Quantum effects can dramatically change the actual 
strength of the link by dressing it with vortex-antivortex 
pairs sitting right at the link] the effective energy of such 
pairs can be quite small [note that a space-time isotropic 
Eq. ([!]) does not apply to RWLs]. However, in accordance 
with the result by Kane and Fisher^ (see also Ref.fTSlfor 
the instanton interpretation) the macroscopic prolifera- 
tion of vortex-antivortex pairs on a link, no matter how 
weak, takes place only at K < 1 (in the thermodynamic 
limit). Hence, as long as K > 1, the quantum physics of 
RWLs can be absorbed into their renormalized Joseph- 
son couplings, a circumstance which we rely upon in our 
analysis, which consequently closely follows ideas devel- 
oped in Ref. [TTJ 

To proceed, we have to specify the distribution of 
RWLs. Here we confine ourselves with a (rather typical) 



strong-disorder setup when RWL emerges as a sequence 
of I* statistically independent microscopic "insulating" 
regions, in which case the probability of finding the RWL 
somewhere in the system of the (dimensionless) size L is 
given by 

Prob[4,L] oc (L/4)e- Cl£ * , (3) 

where e~ Cl defines the probability of a single "insulating" 
region to occur at a given point. We are interested in the 
case when this probability is of order unity, from which 
the dependence of on system size L to logarithmic 
accuracy follows, 

4 ^^(lnL-lnlnL...). (4) 

The second term is related to the Kinchin-Kolmogorov 
law of the iterated logarithm. The Josephson coupling J 
over the link is 

J(4) = F(4)e- C2£ * « F(lnL) (lnL/L) C2 / Cl , (5) 

where the function F(lnL) and the constant e^ absorb 
all the microscopic physics of the weak link. The inverse 
of the superfluid stiffness n s (L) for the system of size 
L ^> £o is then (classically) given by the sum of inverse 
Josephson couplings of RWLs, 

n~\L) = n-\i Q ) + (Jf 1 + J^ 1 + J3- 1 ...)/L. (6) 

As it will be seen self-consistently, the nature of the flow 
at the classical critical point, implied by Eqs. ([5| and (|6|, 
is such that doubling the system size leads only to a minor 
change in n s (L), that is the J value of the weakest link 
gets progressively larger than n s (L)/L as L — > oo. Hence, 
a substantial change of n s (L) happens over a broad range 
of length scales involving a large number of RWLs at 
different length scales. This allows one to treat n s (L) as 
a slow continuous function of L and replace Eqs. ([5|-(|6| 
with the flow equation 

^ = F(ln L)L< + 0(ln In L) , (7) 

where ( = c^jcx — 1. Note that the correction terms 
are iterated logarithms in the system size. By the same 
token, the flow of the compressibility with InL is leaving 
it constant to this order of accuracy. 

The critical condition corresponds to ( —> 0. At 
this point n s (L) = n s (lnL) is decaying logarithmically 
slowly justifying self-consistently the assumption of be- 
ing a smooth function of L. For \(\ <C 1, a substantial 
part of the RG flow remains slow and smooth even on 
the insulating side ( > 0. 

The flow Q predicts K(( = 0, L = oo) = 0. Quan- 
tum effects do not change this flow even for K < K c 
until L is large enough for free instantons to appear. 
The actual quantum critical line corresponds to the 
regime when the classical critical flow takes the system 
to K(L = oo) = K c . If Kq ^> 1, one can reveal how the 
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classical-field flow ultimately crosses over to the univer- 
sal quantum one only at astronomically large L leaving 
no chances to observe it in a cold gas experiment nor 
in first-principles simulations. What can be observed, 
however, is the anomalous finite-size behavior of the su- 
perfluid stiffness. On the superfluid side (£ < 0), finite 
size corrections n s (L) —n s (oo) oc 1/L^I can, in principle, 
be used to extract the critical parameters of the Hamil- 
tonian. 

On the insulating side, a typical feature is the devel- 
opment of substantial sample-to-sample fluctuations of 
n s (L) reflecting the growing effect of single RWLs on 
n s (L). Somewhat away from the critical condition ( = 
(still considering ( < 1) RWLs ultimately destroy the 
self-averaging of n s (L) [and thus applicability of Eq. ^} 
at large enough L. One can imagine that in a peculiar 
regime with condition < £ <C 1 the self-averaging of 
n s (L) survives till K(L = L*) = 1, while the suppressed 
fugacity of instantons prevents the system from having 
insulating properties. In such a case, the Kane-Fisher 
physics on a single RWL will be responsible for render- 
ing the system insulating at L > L*. 

The theory developed in this section can be used for 
processing numeric and/or experimental data, allowing 
one to predict the behavior of the system in question at 
much larger sizes and, in particular, extract critical pa- 
rameters of disorder. A practical example will be con- 
sidered in Sec. |IV[ preceded by Sec. Ill in which we 



address delicate aspects of self- averaging characteristic 
of the strong-disorder classical-field flow, and formulate 
appropriate statistical observables to be associated with 
(marginally) self-averaged quantities of the present sec- 
tion. 



III. STATISTICAL SUBTLETIES. THEOREM 
OF CRITICAL SELF- AVERAGING 

Introduction of the notion of superfluid stiffness as an 
intensive quantity requires self- averaging in the super- 
fluid phase, at least in a certain minimalistic sense re- 
quiring that relative fluctuations of n s (L) among differ- 
ent disorder realizations be small and get progressively 
smaller with increasing L (see also further on in this sec- 
tion and in Sec lIVI for a numerical estimate of the mini- 
mum system size required for this to hold in the super- 
fluid phase). The same is true for n^T^L), the prime ob- 
ject of the analysis of the previous section. Nevertheless, 
at small the issue of self- averaging becomes delicate 
because of the potential divergence of the dispersion of 
nj 1 ^) in the L — >• oo limitPM [This divergence does 
not necessarily contradict with a vanishing characteristic 
width of the distribution.] 

While not creating fundamental theoretical difficulties, 
the divergence of dispersion requires some care in pro- 
cessing finite-size data. Obviously, one has to avoid work- 
ing with the dispersion when characterizing fluctuations 
of n^^L). Correspondingly, in order to characterize the 



distribution of a random variable x we use three typical 
numbers (below P stands for probability): (i) the median 
x such that P[x < x] = 1/2, (ii) the lower characteristic 
value X- such that P[x < X-\ — 1/4, and (hi) the upper 
characteristic value x+ such that P[x > x+] = 1/4. The 
width of the distribution is then naturally characterized 
as Sx = x+ — X-. Correspondingly, by the relative width 
we understand R = Sx/x. 

A special convenience of working with x± and x is due 
to the fact that for any random number y = f(x) related 
to x by a monotonic function /, we have y = f(x) and 
V± = f( x ±): or U± = I( x t)' I n wna t follows all physical 
quantities are characterized by their medians and (rela- 
tive) widths, not the first and the second moments of the 
distribution functions. 

We now prove a theorem that self-averaging in the 
above mentioned sense takes place at the critical line 
separating SF and BG phases, no matter whether it is of 
the GS type, or corresponds to some hypothetical strong- 
coupling scenario based on weak links, which features a 
critical value of the Luttinger liquid parameter K c > 3/2. 
Finite K c implies finite critical value of 1 /n s because with 
the compressibility remaining finite at the critical point 
we can speak of K and l/n s interchangeably at the level 
of their medians. This is consistent with the most conser- 
vative definition of the superfluid phase and the critical 
point at its boundary as states with non-zero probability 
density, P^ L \x), of finding a finite value of x = l/n s in 
the thermodynamic limit L — >> oo. This leads to a mean- 
ingful (i.e. with non-zero weight over a large-enough fi- 
nite interval) limiting function which is allowed — but not 
a priori required — to be a ^-function: 



P{x) = lim pM(x). 

L— ^oo 



(8) 



Consider now the process of combining two macroscop- 
ically large systems of size L into one system of size 2L 
by adding a link between their edges. Due to the local 
renormalization of the weak link physics we accept that 
the contribution of the link to the macroscopic stiffness 
may depend on properties of the left and right pieces, 
but not on their macroscopic size. In other words, if on 
the left the superfluid stiffness is x and on the right it 
is then the link distribution function is Q x , y (ti = 1/J) 
and the composition law for the probability distribution 
of z = l/n s in the combined system is [see ([6|] 



pOO /'OC r& 

p( 2L )(z) = dx dy 

Jo Jo Jo 



x+y+Z/L 



(9) 



The distribution of weak links is certainly normalizable. 
Moreover, even at classical criticality it decays at large 
£ as a Lorentzian. Since, at any fixed z, the ranges of 
integrations in (I9J) are finite, we can safely take the limit 
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of L ^ oo and replace P with P 



P(z) 



dx dy dt; 
o Jo Jo 



P(x)P(y)Q X)y (OS 



x+y+t/L 



Integration over z in Eq. ( 10 ) then yields 
f = fq, 

where 



(10) 



(11) 



/ 



PCX) 

= / P(x) 

JO 



dx < 1 . 



/ 0(0 

Jo 



df=l. (12) 



Since / > (without a finite fraction of systems with 
non-zero n s , we fail to meet the definition of the critical 



state with finite n s ), from Eq. (11) we immediately con- 
clude that / = 1, meaning that at the critical point the 
limiting distribution P(x) is normalized to unity. 
In Fourier space, Eq. (fTob reads 



d# dy P(x)P(y)e 



ik(x+y)/2 



I 



ifc£/2L 



(13) 

Since the last integral corresponds to macroscopically 
small Fourier harmonics of Q XjV that decay at large £ 
at least as fast as a Lorentzian, we can replace it with 
unity [the leading correction is 0(1/ L)]. This leaves us 
with the equation 
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(14) 



for which we have two options: (i) Either = 1, or 
(ii) |P/c| decays exponentially at \k\ —> oo. Option (ii) 
is incompatible with the requirement that the function 
P(x) be identically zero at x < 0. Hence, we have = 
e %k/n s ^ • e p(#) = S(x—l/n 3 ), which proves the theorem. 

Two remarks are in order to prevent possible confu- 
sions. Having observed that typical values of J in the 
composition law ([9| are of order unity, while the strength 
of characteristic weak links responsible for, say, the clas- 
sical flow Q is macroscopically small, one may wonder 
how the two facts should be reconciled with each other, 
and, in particular, how those macroscopically small links 
are properly taken into account in the composition law 
([9]). The first confusion is removed by noting that the 
weak links in Eq. ^ are being selected throughout the 
whole system, while the link in the composition law ([9| 
is found at a given position in space, right at the center 
of the system. The macroscopic difference in two pro- 
cedures thus accounts for the macroscopic difference in 
characteristic values of the corresponding links. The role 
of macroscopically weak links is implicit in Eq. (|9| — 
those links effect the form of the distribution P^Jx). 

Our second remark concerns confusing self- averaging 
with the existence of a finite dispersion. One should not 
think that P(x) = S(x — l/n s ) somehow implies a finite 



dispersion at any finite L. The dispersion depends on 
the limit taken to obtain the generalized function. Self- 
similar distributions for random variables based on sums 
of identical random v ariabl es are known in mathematics 
as stable distributions! 22 1 23 * A distribution for a random 
variable X is said to be stable if for two independent 
copies Xi and X2 and for any constants a > 0, b > 
it holds that aX\ + 6X2 = 7X + d for some constants 
7 > and d. Normal distributions are stable distribu- 
tions (by the central limit theorem, when the variance is 
finite), but there exists a generalized central limit theo- 
rem which states that the sum of a number of random 
variables with power-law tail distributions decreasing as 
where < a < 2 (and therefore having infinite 
variance) will tend to a stable distribution of the from 
P(x) = ^ J e ltx e~( 7t ^ dt, with scale parameter 7 > 
which is a measure for the width. Note that this goes 
beyond the law of the large numbers. The distribution 
can be centered around \i by the substitution x —> x — /i. 
In particular, for a = 1 we have the well-known Lorentz 
distribution which has infinite moments at any finite 7. 
Physically, we know that the divergence of the disper- 
sion with L is certainly happening for a marginal super- 
fluid, but the system will ultimately flow to an insulator. 
However, one can also envisage situations in a superfluid 
where the bi-modality due to the divergence of the dis- 
persion is resolved by keeping the location of the center 
Ii finite while the statistical importance of the tail is ir- 
relevant because of macroscopical pre-factors suppressing 
the tail (cf. the stable distributions above). Such exam- 
ples should not be thought of as a different superfluid 
phase (from the usual one with vanishing dispersion at 
any finite L) because self-averaging does take place in the 
thermodynamic limit with only rare, thermodynamically 
unimportant fluctuations. 



The theorem we proved can, in particular, be used to 
identify which systems flow to insulators (cf. Fig [4|. 
Moreover, self- averaging on the critical line allows one 
to use the hydrodynamic description at large distances. 
By continuity, even if one increases disorder strength be- 
yond the critical value by some arbitrary small but finite 
amount, A = A c + e, the hydrodynamic description will 
apply up to enormously large distance scales. Next we 
notice (by the instanton analysis described above), that 
if this is taking place with K c > K > 3/2, then disorder 
fluctuations become less and less relevant with increasing 
the distance scale, i.e., the state will retain its superfluid 
properties and will not become an insulator, in formal 
contradiction with the K c > 3/2 scenario. We do not 
see any physical mechanism which invalidates this con- 
clusion. 
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FIG. 3. (Color online). Median inverse stiffness, n^ 1 , and 
compressibility, k, as functions of system size L for disorder 
strength A = 1.4. (Compressibility error bars are shown but 
too small to be resolved in this plot.) 



IV. ILLUSTRATIVE SIMULATION 

A. Model and Technique 

We consider the disordered J-current modeP^, 
^2tJi,x ~ \U\Jl iT ~ (/x - e x )J niT . 



100 



1000 



H 



(15) 



Here, the integer vector n = (x, r) labels sites of the 
two-dimensional square space-time lattice of dimension 
L x x L r , J n ^ a are integer "currents" living on lattice 
bonds labeled by n, a where a = x,r stands for unit 
vectors pointing in the positive discrete space and imag- 
inary time directions, respectively. The allowed values 



for the currents are \J n 



< 1 and < J n , T < 2. 



(We 



will write L x and L intermittently for the system size). 
The allowed configurations of bond currents are subject 
to the zero-divergence constraint ^ a (J n ,-a + 4,a) = 0? 
where, by definition, J n -a = —Jn-a,a- We consider at- 
tractive interactions and choose them as the unit, U = 1. 
Disorder e x uniformly chosen in the box [—A, A] is added 
to the non-zero chemical potential, which breaks particle- 
hole symmetry (the SF-BG transition is then expected at 
K c = 3/2 according to the GS scenario). In the absence 
of disorder, a first order SF-MI transition occurs in a 
broad parameter range starting from (fi/U = — 2,£ = 0). 
To test the theory on presently achievable system sizes, 
we take t = 1/4 and \i = —1.8 in the vicinity of this 
first-order transition, and expect strong-disorder fluctu- 
ations to generate local insulating regions responsible for 
RWLs. 

The central (dimensionless) quantity of interest is the 
superfluid stiffness, n s . Also important is the compress- 



FIG. 4. (Color online). The relative width R s = n s Sn^ 1 of 
the distribut ion of n^ 1 (calculated in accordance with conven- 
tions of Sec. |HI| ). The increase in relative width for A > 1.28 
leads to insulating behavior, according to the theorem of self- 
averaging on the critical line. 



ibility, ft, since a combination of n and n s gives the Lut- 
tinger parameter 

K = it^n s K . (16) 

The superfluid stiffness and compressibility can be com- 
puted from the winding numbers W x and W T in space 
and time, respectively: 

n. = «^ 2 > " (W x f) (L X /L T ) , 



K =((W?)-(W T ) 2 )(L T /L X ). 



(17) 
(18) 



Here, the angular brackets stand for statistical averaging 
for a given disorder realization. Assuming a transition 
with dynamical exponent z = 1, we perform finite size 
scaling with constant aspect ratio L T /L X = 1/2. In the 
figures we plot results for winding number fluctuations 
which explains an extra factor of 1/2. Simulations were 
done by employing the classical worm algorithm 25 . We 
typically consider between 5,000 and 20,000 disorder re- 
alizations in order to have converged answers; from a 
couple of 100 realizations one cannot determine the dis- 
tribution functions. 



B. The Results 

Our first observation is that, in accordance with the 
theory of Sec. |n| (as well as with the standard GS sce- 
nario), it is only the superfluid stiffness that gets sub- 
stantially renormalized in the critical region and in the 
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FIG. 5. (Color online). The flows of the median values of 
n~ x (L) fitted with Eq. |7| in which the (unknown) function 
F(lnL) is chosen to be a (A-dependent) constant, F = B(A), 
treated as a fitting parameter (see also Fig. [7|. 



FIG. 6. The exponent C(A). The arrow indicates the clas- 
sical critical value of disorder, A^ = 1.275(10), such that 
= 0. The function C(A) is structureless in the vicin- 
ity of A = A^ cl) . 



insulating phase, while the compressibility remains con- 
stant (within error bars) for all system sizes, see Fig. [3] 

In the presence of strong disorder^the distribution of 
n^ 1 becomes quite broad (see App. [A|for the full distri- 
butions) . The flows of relative widths as functions of the 
system size, R S (L), are shown in Fig. H In the theoret- 
ical limit of large L, the function R S (L) is supposed to 
decrease (tending to zero at L — >• oo) in the superfluid 
phase, and increase in the insulating phase. As we see in 
Fig. |4j these two distinctively different asymptotic trends 
set in only at L > 100 (with the separatrix at A < 1.28). 
At L < 100, the relative width R S (L) increases on both 
sides of the separatrix. This, however, does not mean 
that the theory of Sec. [il] does not apply at L < 100. 
Actually, the opposite is true since the initial growth of 
R S (L) is naturally attributed to the fact that the func- 
tion F(lnL) in Eq. Q is small already at L ~ 1, rather 
than only in the asymptotic limit of L ^> 1. Indeed, 
as is clear then from Eq. Q [that now applies almost 
from L ~ 1] 5 for n" 1 to get significantly renormalized — 
which is a necessary condition for developing substantial 
sample-to-sample fluctuations — the system size L has to 
increase dramatically. 

Figure [5] reveals a remarkable quantitative agreement 
between the numerical results for the median values of 
n7 x (£) and the theory of Sec. [nj Our fitting procedure 
is based on Eq. Q in which the function F(ln L) is chosen 
to be a (A-dependent) constant, F = 5(A), treated as 
a fitting parameter for a given strength of disorder. The 
second fitting parameter is the exponent C(A). A crucial 



consistency condition then is that the functions B(A) 
and C(A) be structureless in the vicinity of the classical 
critical value of disorder, meaning that our data reveals 
nothing beyond Eqs. (|3|-([7|). As we see in Figs. [6]and[7| 
this condition is met. 

The genuine (quantum) critical value of disorder, A c = 
1.265(10), is found by extrapolating K(A,L) to L — 
oo and requiring if(A c ,oo) = K c = 3/2, see Fig. [8] 
This procedure completely ignores the renormalization of 
K(L) due to instanton pairs, which is justified by the van- 
ishingly small fugacity of the latter, as argued in Sec. [TTJ 
The value of A c has overlapping error bars with the clas- 
sical critical value A^ = 1.275(10) obtained in Fig. ^ 
This is due to the fact that K is a steep function of n s 
at n s — > 0, which is the reason we can find extrapolated 
if (oo) values very close to K c = 3/2. [Our error bars 
do not account for possible systematic bias originating 
from the simple choice of F = B(A). The corresponding 
analysis, though possible, goes beyond the goals of the 
present paper.] 

Figure |9j to be compared to Fig. [T] and contrasted to 
Fig. [2] shows the flows of K(L). The qualitative agree- 
ment with Fig. [l] and the pronounced difference from 
Fig. [2] are clearly seen. Especially instructive are the 
data sets for A = 1.6 and A = 1.75. At those val- 
ues of disorder the system is unquestionably well inside 
the insulating phase — to the extent that even fitting the 
flow by the ansatz of Eq. ([7]) can hardly be justified be- 
cause Eq. ([7]) is derived under the assumption of ( <C 1, 
while C(A = 1.6) w 0.8 (see Fig. [6]). Nevertheless, the 



8 



2B(A) 



I 



0.2 



5 



0.1 



0.9 



1.1 1.2 1.3 1.4 1.5 



K(oo) 

3 
2.5 



1.5 



1 



0.5 




0.95 1 1.05 1.1 1.15 1.2 1.25 1.3 A 



FIG. 7. The parameter F = 5(A) extracted from the fits. 
Within the error bars, the dependence of B on A is linear, as 
shown by the dashed line. 



marginal superfluid behavior persists at A = 1.6 and 
A = 1.75 till system sizes of few hundreds. Most impor- 
tantly, the flows for A = 1.6 and A = 1.75 go well below 
the quantum-critical value K c = 3/2 without any change 
of their character. Qualitatively, the flows for A = 1.6 
and A = 1.75 are very similar to the flow for A = 1.5 ex- 
trapolated to K(L) ~ 1 5 which is dramatically different 
from the situation sketched in Fig. |2j 

Note that for A = 1.6 and A = 1.75 marginal super- 
fluidity persists till the Kane-Fisher point Kkf = 1, in 
agreement with our conclusion concerning the suppressed 
fugacity of instanton-anti-instanton pairs. Because for 
our model the classical separatrix and quantum critical 
line are so close to each other, there is no chance of see- 
ing a quasi-plateau between K = 1 and K = 1.5. The 
slow flow thus shows that, on physical length scales, we 
are dealing with a mesoscopic situation where the notion 
of the thermodynamic state is less meaningful than the 
flow itself: Making cuts at fixed L or fixed K = 2.5 in 
Fig. [9] re veals little information in itself but is meaningful 
when understood in terms of the flow equations. Hence 
accurately tracing the evolution of K with L is abso- 
lutely crucial for understanding the underlying physics. 
A recent numerical study 26 of a particle-hole symmetric 
(l+l)-dimensional XY model came to opposite conclu- 
sions as the ones presented here. The authors of Ref. l26l 
analyzed their data by monitoring the susceptibility, Lut- 
tinger parameter, and correlation length £ and found a 
non-universal K c for strong enough disorder. Their type 
of analysis is however incompatible with our data, but 
we believe that resolving the flow of K with L on the 
critical separatrix for strong disorder, which is lacking in 



FIG. 8. The asymptotic values of the Luttinger parameter 
extracted from extrapolations of classical flows to L — >> 00. 
The dashed line is the fit K(A,L — 00) = ay/x/(x + 6), 
x = A ( c cl) -A (a and b are fitting parameters) consistent 
with the behavior of K(A,L = 00) predicted by Eq. |7]). 
The quantum-critical value of disorder, A c = 1.265(10), is 
obtained from the condition K(A C , L = 00) = K c = 3/2. 



the paper, could settle the discrepancy in favor of the 
theory outlined in this work. 



V. CONCLUSION 

Our work proves that strong disorder may lead to a 
logarithmic flow of the stiffness coming from the classical- 
field renormalization while compressibility remains con- 
stant. Although the universality class of the superfluid 
to Bose glass quantum phase transition does not change, 
the classical-field mesoscopic physics renders this univer- 
sality essentially non-observable on physically relevant 
length scales, since a peculiar mesoscopic behavior per- 
sists up to extremely large system sizes. Even deep in 
the insulating side a physical (e.g., cold-atomic) system 
can demonstrate a marginal superfluid response, with 
strong sample-to-sample fluctuations. It makes hence lit- 
tle sense to identify the thermodynamic phase for certain 
Hamiltonian parameters if only finite samples are avail- 
able. We believe that this mesoscopic behavior, illus- 
trated by numerical data, is the long-sought reconcilia- 
tion of classical-field criticality with the established pic- 
ture of the superfluid-to-glass quantum phase transitions 
(which remains intact in the academic thermodynamic 
limit). Fitting numerical/experimental data for the su- 
perfluid stiffness as a function of disorder and the system 
size with the classical-field renormalization flow allows 
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FIG. 10. (Color online). Distribution of the compressibility 
k for A = 1.28 close to the critical line. 



FIG. 9. (Color online). The flows of median values of the 
Luttinger parameter. In view of a perfect self- averaging (and 
absence of any appreciable L-dependence) of k, the solid lines 
correspond directly to their counterparts in Fig. [5] processed 
in accordance with Eq. ( |16[ ). The dashed lines for the A = 1.6 
and A = 1.75 data are added simply to guide an eye (see the 
discussion in the text). 

one to accurately extract the critical value of disorder 
and predict the behavior of the system at much larger 
sizes. In future work, we plan to employ this theory for 
completing the SF-BG phase diagrams of the disordered 
Bose-Hubbard model in the weakly interacting Bose gas 
regime in three^, twcP, and one dimension. 

Finally, we proved that — regardless of the possible di- 
vergence of variance of inverse superfluid density — any 
quantum critical flow possesses self- averaging properties 
and is thus subject to an asymptotic hydrodynamic de- 
scription. In our view, this theorem rules out quantum 
strong-coupling criticality scenarios. 

We are grateful to Ehud Altman, Thierry Giamarchi, 
Victor Gurarie, David Pekker, Susanne Pielawa, and 
Anatoli Polkovnikov for valuable discussions. This work 
was supported by the National Science Foundation under 
the grant PHY-1005543. LP thanks the hospitality of 
the Kavli Institute for Theoretical Physics at UCSB and 
at Beijing, where part of this work was done. 



Appendix A: Full Distribution Functions 

Although it is argued in this work that the median and 
the (relative) width defined in terms of percentiles of the 
distribution give a perfect description of the physics, we 
nevertheless want to show the full distribution functions 



0.12 




FIG. 11. (Color online). Distribution of the stiffness n s for 
A = 1 in the superfluid phase. Arrows indicate the 25, 50 
and 75 percentile, from left to right for L — 16, L — 64 and 
L = 256. 



in this Appendix. On the one hand this gives the reader 
a better appreciation of the broadness of these distri- 
butions (and the corresponding numerical challenge to 
resolve the distribution) as well as an indication of where 
the median lies. On the other hand, the strong disorder 
scenario of Ref. [18] puts special emphasis on the tails of 
the distributions for the stiffness n s and we also would 
like to provide this information for completeness. 

We start with the distribution for the compressibil- 
ity. As argued in the text, the compressibility plays no 
important role in the flow. In Fig.[l0]we show the distri- 
bution of k for different system sizes at A = 1.28 which 
is very close to the critical line. The distributions show 
strong self-averaging when increasing the system size, as 
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FIG. 12. (Color online). Distribution of the stiffness n s for 
A = 1.28 close to the critical line. Arrows indicate the 25, 50 
and 75 percentile, from left to right for L = 16, L = 64 and 
L = 256. 



broad for small system sizes). 

We proceed with the distributions of the stiffness n s . 
Deep in the superfluid phase we see with the naked eye in 
Fig. 11 for A = 1 that self- averaging occurs even though 
it is a slow process with system size. The median tracks 
the maximum of the distribution well (they will coincide 
in the thermodynamic limit). When we increase A to be 
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FIG. 13. (Color online). Distribution of the stiffness n s for 
A = 1.4 close to the critical line, but in the insulating regime. 
Arrows indicate the 25, 50 and 75 percentile, from left to right 
for L = 16, L = 64 and L 256. 



expected. So for k we are able to define the average 
and the variance without problem. This picture is not 
altered in the superfluid or in the insulating phase (not 
shown). We are thus safe to ignore the renormalization of 
the compressibility (even though the distribution is very 



FIG. 14. (Color online). Distribution of the stiffness n s for 
A = 1.6 in the insulating regime. Arrows indicate the 25, 50 
and 75 percentile, from left to right for L = 16, L — 64 and 
L 256. 



close to the critical point (see Fig. 12 for A = 1.28) the 



median shifts to lower values for small system sizes. How- 
ever, the flow seems to converge to a finite value in the 
superfluid phase, while the relative width is decreasing 
in the superfluid phase and increasing in the insulating 
phase (see Fig. [4|. On the numerically accessible sys- 
tem sizes, nothing is seen to occur around n s = near 
the critical line. For slightly larger A (see Fig. 13 for 
A = 1.4) the flow to the insulator can nicely be tracked. 
For small system sizes such as L = 16 the distribution 
has a well- developed peak but we note an asymmetric 
and pronounced tail for smaller n s that develops strongly 
when we increase L. For L = 256 the distribution ap- 
proaches n s = in such a way that the variance for nj 1 
certainly diverges. After that, the insulating behavior is 
unstoppable and quickly becomes apparent in a similar 
fashion as seen deep in the insulating regime (see Fig. 14 
for A = 1.6) already at smaller system sizes (L = 128 
and L = 256): The distribution quickly broadens and 
develops a strong overlap with n s = for increasing L. 
Ultimately the median will flow to 0. 
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